Constraining the Variation of G 
by Cosmic Microwave Background Anisotropies 
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We use the Cosmic Microwave Background Anisotropies (CMBA) power spectra to constrain the 
cosmological variation of gravitational constant G. It is found that the sensitivity of CMBA to the 
variation of G is enhanced when G is required to converge to its present value. The variations of G 
from the CMB decoupling epoch z ~ 1000 to the present time are modelled by a step function and 
a linear function of scale factor a respectively, and the corresponding 95% confidence intervals for 
G/Go are [0.95, 1.05] and [0.89, 1.13], Go being the present value. The CMBA constraint is unique 
in the sense that it entails the range of redshift from z « 1000 to 0. 
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I. INTRODUCTION 

One of the fundamental questions in physics is whether 
the fundamental constants are truly constant. Indeed 
the possibility of cosmological variation of "constants" 
has long been proposed d]. Among the fundamental con- 
stants, the gravitational constant G is the least accu- 
rately measured. The value of G is measured in the 
laboratory and applied to all scales. To check for the 
constancy of G, tests should be done at different spatial 
and temporal scales. There are many tests at redshifts of 
order 0. For example, the lunar laser ranging experiment, 
which monitors the distance between the earth and the 
moon by laser ranging technique, can be used to put a 
bound on the variation of G0] ■ If G varies during the his- 
tory of the earth, the surface temperature and size of the 
earth would change. However the earth does not preserve 
a good record of the gravitational conditions. Increase in 
G causes the Sun to burn at a faster rate, and the depth 
of the convective zone is affected, which can be observed 
in the vibrational modes of the Sun, in particular the p 
mode 3. Since an increase in G also shortens the life 
spans of stars, the ages of stars in globular clusters can 
be used to put a constraint on the deviation of G from 
its present value M. Because the Chandrasekhar mass 
Mch oc G~^/^m^ sets the mass scale of neutron stars, 
by observing the masses of neutron stars formed at dif- 
ferent redshifts, limits on G in the past 10 Gyr or so have 
been derived in [5]. The highest redshift, of about 10^", 
constraint comes from Big Bang Nucleosynthesis (BBN) . 
An increase in the expansion rate during the epoch of 
BBN causes the freeze-out to occur earlier, and the abun- 
dances of neutron and hence ''He are enriched @. For 
details of various experiments and observations, see the 
review article by UzanQ (see also [1]). 

On the theory side, there have been grand unification 



theories and string/M theory motivated models predict- 
ing that some of the fundamental "constants", such as 
the fine structure constant a and the Newtonian gravita- 
tional constant G, may vary over time. In theories with 
extra dimensions, the effective gravitational constant in 
4D spacetime depends on the more fundamental mass 
scale in the bulk and the size of the extra dimensions^. 
If the size of the extra dimensions evolves over time, the 
effective constants in 4D will vary as well. For example, 
the DGP model 10] has been put forward to explain the 
recent observation that the universe is in an accelerating 
phase without invoking the dark energy. It was argued 
that the acceleration is due to the leakage of gravity into 
the extra dimensions. 

To select promising ones from the myriads of models 
in the literature, one may constrain possible variations 
of the fundamental "constants" using observational data. 
There have been efforts trying to constrain the possible 
variation of G using cosmological data. In particular, 
since the Cosmic Microwave Background Anisotropies 
(CMBA) is sensitive to many cosmological parameters, 
it could be used to constrain the variation of G. CMBA 
is unique because it offers a long look back time. The 
physics of CMBA is particularly clean as it involves only 
well known physics in the linear regime. One approach is 
to constrain the variation of G in some particular types 
of models, e.g. the CMBA spectra in the Brans-Dicke 
cosmology are discussed in Ref. [HI, HI]- However, given 
the multiplicity of models in the literature, it seems more 
practical to use a simple and generic parametrization for 
G. In [H, [11], the authors have used the CMBA power 
spectra to constrain the possible variation of G with a 
parameter A as 
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where Go is the present laboratory-measured value. It 
was assumed that A is a constant over the age of the uni- 
verse (and only suddenly becomes 1 at the present time) . 
However this assumption is unrealistic since we know that 
G should converge to its present value to avoid confiicts 
with other experiments and observations. We shall call 
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the convergence of G to its present value stabilization. 
One can imagine that G may vary in many different man- 
ners over the history of the universe, and so it is hopeless 
to deal with all possibilities one by one. In this article 
we study two generic stabilization schemes. One of them 
is instantaneous stabilization. That is we consider an 
abrupt gravitational transition and model the variation 
of G by a step function: 

^2^rA2ifa<a„ 

[ 1 li a > fls , ^ ' 

where a is the scale factor and is the scale factor 
at which stabilization occurs. Another is that G varies 
smoothly and we parametrize it as a linear function of a: 

{Ag if a < a*, 

l-^(l-^o) ifa*<a<«. , (3) 
1 if a > Os , 

where is the scale factor at the time of photon decou- 
pling. Our main goal is to constrain the range of Aq in 
Eq. [5] and Eq. [31 We assume that the underlying mecha- 
nism for the variation of G does not affect other physics 
so that the standard CMBA calculation with the modi- 
fications of Eq. [1]- [3] is valid. The Friedmann equation is 
modified as 

{^^^' =X\a)H^{a), (4) 

with 

H\a)^Hl(^ + ^ + n^, (5) 

where a dot denotes derivative with respect to the con- 
formal time, Hq is the present Hubble parameter, U,m is 
the density parameter of the non-relativistic matter, f2^ 
is the density parameter of radiation, and I^a is the den- 
sity parameter of the cosmological constant. Note that 
we consider a flat universe here. 

The paper is organized as follows. The effects of vari- 
ation of G on the CMBA power spectra are investigated 
in Section ini In Section Hill we study the constraints on 
the variation of G by the three year WMAP data using 
the method of Markov Chain Monte Carlo (MCMC) and 
discuss the results obtained. Section HVl is devoted to the 
conclusion. 



II. EFFECT OF VARIATION OF G ON THE 
CMBA TEMPERATURE AND POLARIZATION 
POWER SPECTRA 

It has been pointed out in Ref. [l^ that the CMBA 
angular power spectrum does not change using the sim- 
ple prescription Eq. [H as far as the dynamical equations 
are concerned. The ratio between the sound horizon and 
the distance to the last scattering surface (LSS), 8, will 
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FIG. 1: The CMBA temperature angular power spectra for 
three values of the Newtonian gravitational constant with 
instantaneous stabilization Zs = 0. The solid, dotted and 
dashed curves correspond to Aq = 1, 0.5 and 3 respectively. 
The damping effect is only noticeable for I > 500. 

not change if both are blown up by the same factor due 
to the varied gravitational constant in the flat universe. 
However, the scaling is not perfect since the recombina- 
tion physics does introduce another scale. The recombi- 
nation is dictated mainly by the binding energy of hy- 
drogen atom, which is not affected by a change in the 
gravitational constant. If the expansion rate is greater 
during the epoch of recombination, it will be more diffi- 
cult for the protons and electrons to recombine, and the 
ionization fraction Xe will increase. The probability den- 
sity that a photon last scatters at a conformal time rj is 
given by the visibility function (for a review of CMBA 
physics, see for example [l6|), 

g[n) = fe~^, (6) 

with 

f = aUeCTT, (7) 

where Ue is the number density of electrons and ctt is 
the Thomson cross-section. The increase in Xg broad- 
ens resulting in more severe damping of the high 
I peaks. However the duration that the photons stay in 
contact with the LSS is also shortened as the expansion 
rate is greater. Because the two effects partially can- 
cel each other, the damping is not increased much even 
when G is increased by a large amount. In Fig. [1] we 
show the temperature power spectra for Aq = 0.5 and 3 
compared to the unchanged one (Ag — 1) for the case of 
of instantaneous stabilization with the stabilization red- 
shift Zs = — — 1 = 0. A large change in G is required for 
noticeable changes in the spectra. The damping effect 
can be partially compensated by reionization, and so the 
details of reionization such as the degree of reionization 
affect the sensitivity of the spectra to A. 

The above conclusion that the CMBA angular power 
spectrum is not sensitive to the value of G is based on 
the assumption that the gravitational constant remains 
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FIG. 2: Same as Fig. [T] but with spectra with Zs = 10 in 
contrast with those with Zs = 0. The standard spectrum 
(solid) nearly coincides with the one with Aq — 1.2 and Zs ~ 
(dot-dashed). When instantaneous stabilization occurs at 
Za — 10, marked shifts to larger (smaller) I scales result for 
Ag = 1.2 (dotted) (0.8 (dashed)). 



FIG. 3: The E-polarization spectra with instantaneous stabi- 
lization at = and 10. The peaks shift to larger (smaller) I 
scales for Aq = 2 (0.6) when instantaneous stabilization takes 
place at Zs — 10, in contrast to the cases without stabilization 
{Zs = 0). 



different from the present laboratory-measured-value till 
"yesterday." However G must converge to its present 
value so that there is no conflict with the low-redshift 
constraints. When G changes as described in Eq. [2] or 
Eq. [21 the expansion history of the universe is modified, 
the fractional change in the sound horizon at the epoch 
of decoupling is different from that in the conformal dis- 
tance to the LSS, and the resultant CMBA spectra will 
be distorted. Fig. [5] shows the temperature power spec- 
tra with instantaneous stabilization at = and 10. 
Note that the "standard" spectrum (Aq = 1 and Zg = 0) 
nearly coincides with the one for Aq — 1.2 and Zg = 0- If, 
however, the stabilization redshift is at Zs = 10, the spec- 
trum for Ao=1.2 (0.8) shifts to larger (smaller) / scales. 
In fact, it is conceptually simpler to compare the one for 
A§ = 1.2, = to the one with Aq = 1.2 and Zs — 10. 
The size of the sound horizon is the same for both cases 
while the distance to the LSS is increased for Zg — 10; 
as a result 8 becomes smaller and the spectrum shifts to 
high I scales. The same argument applies to the one with 
Ag = 0.8 and Zg = 10, but with the opposite effect. Be- 
cause of the dramatic gravitational transition, we observe 
an enhanced late Integrated Sachs Wolfe (ISW) effect in 
the small I scales. 

We observe similar sideway shifts in the E-type po- 
larization and TE cross polarization spectra as well. In 
Fig. [3l we show the E-polarization power spectra with in- 
stantaneous stabilization at Zs = and 10. In Ref. p^ . 
the authors proposed that the degeneracy between the 
expansion rate and the scalar spectral index Us can be 
lifted by measuring the polarization, because the forma- 
tion of polarization is proportional to the width of the 
visibility function. An increase in G causes the power 
of the polarization spectrum to increase in small I scales 
and then decrease in large I scales. But the effect of 
stabilization is much more appreciable. 

When G varies linearly with a, its effects on the CMBA 
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FIG. 4: The effects of linear variation of G on the temperature 
power spectrum. The solid line shows the spectrum when 
there is no variation of G. When Aq = 1.2 (0.8), the peaks 
shift to larger (smaller) I scales, as is set to be 1. 

power spectra are still dominated by the sideway shifts 
discussed in the instantaneous stabilization scenario, as 
can be seen in Fig. [J] The effects are qualitatively similar 
for the two stabilization schemes that we use. 

Since our prescription is simple, we can easily estimate 
the amount of shift in the CMBA spectrum due to stabi- 
lization. This calculation is similar to that of the CMB 
shift parameter [l5|. From Eq. [4l we have the conformal 
distance to the LSS 

where 770 is the present conformal time, and 77* is the con- 
formal time at CMB decoupling. One immediately sees 
that if Aq is larger, the conformal distance will be smaller. 
This is reasonable because it takes less time for the uni- 
verse to expand to its present size. The peak positions 
in the CMBA temperature angular power spectrum in a 
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flat universe can be characterized by [li 



In « rnr — , 
where rc is the sound horizon: 



61n 



51, 



sdr]. 



The sound speed Cs is given by 

1 



■i{l + R)'' 



with 



4 



(9) 



(10) 



(11) 



(12) 



where pB and p-y are mass densities of baryons and ra- 
diations respectively. It should be pointed out that the 
sound horizon also depends on A through r/. Hence if 
Ao > 1 at the epoch of CMB decoupling, the size of the 
sound horizon will also be smaller. In fact, the reduction 
rates for and are the same if there is no stabilization; 
the effects of A are exactly cancelled, a manifestation that 
the peaks do not shift if there is no stabilization. How- 
ever, when there is stabilization, two different scales will 
be introduced and the peak positions will shift. 

We now illustrate with the case of instantaneous sta- 
bilization. The calculations can be simplified by noting 
that if there is no stabilization, the peaks do not shift 
even if Aq ^ 1. The sound horizon is the same irre- 
spective of stabilization if it takes place after decoupling. 
Hence we only need to compare the conformal distance 
between the case with stabilization and the one without 
it: 



— a*s 
1 



da 



am{a)' 



(13) 
(14) 



where the subscript jvs denotes no stabilization and 5 
denotes with stabilization. From Eq. [9l the shifts in the 
peak positions due to the change in G are given by 



Sin — — Sd^. 



(15) 



Plugging in the standard ACDM model parameters, we 
have 

Sin 5087i(Ao - 1)(0.558- VO.Slla^). (16) 

One may proceed similarly for linear variation of G, but 
it is too cumbersome to write down the results explicitly. 
We display 6ln/n = n6{d^,/rs) against A§ in Fig. [5l We 
see that the shift in I from the simple arguments here 
agrees with the full numerical calculations in Fig. [2] and 
Fig. m Furthermore, the curve by the simplified formula 
in Eq. (TH] tracks closely the one from complete calcula- 
tions of r, and d«. 
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FIG. 5: The left (right) panel shows 5ln/n = Tr5{dt/ra) vs. Aq 
for instantaneous (linear) stabilization. For the plot on the 
left, the thick (red) solid curve corresponds to as=0.1 and 
the dashed curve to as— 0.5. The dots correspond to Sln/n 
averaged over the first three peaks of the spectra in Fig. [2] 
(Ao = 1.2, Zs = 10 and Aq = 0.8, = 10). The two thin 
solid (black) curves are obtained using the simplified formula 
Eq. [16] for = 1 and 0.5 respectively. On the right, the 
solid and dashed curves correspond to linear stabilization with 
as = 1 and 0.5 respectively. Similarly we denote by the dots 
the average 5ln/n in Fig. 2] for as — 1, and Ao = 1.2 and 0.8 
respectively. 



III. CONSTRAINING G BY MCMC AND 
DISCUSSIONS 



Because the CMBA angular power spectra are sensitive 
to various parameters, to be consistent, other relevant 
parameters should also be varied when fitting with data. 
A popular means is to make use of the Markov Chain 
Monte Carlo (MCMC) method. With the temperature 
and polarization spectra computed by the Boltzmann 
code CMBFASTM, we employ the public MCMC en- 
gine CosmoMCpTj to explore the parameter space. Since 
the Hubble parameter Hq is measured to rather good pre- 
cision by the HST key project, we take Hq = 72|JJj]. The 
constraints we get are not sensitive to this restriction. 
Thus the free parameters that we vary are: ujb = ^sh^ , 
i^CDM = ^CDAih^, Zre, the reionization redshift, Ug, the 
index of the primordial perturbation spectrum. As, the 
normalization amplitude, A§, and the stabilization red- 
shift Zg. We use the three year WMAP datajl^l to con- 
strain these parameters. 

First of all, we do not consider stabilization; that is, we 
assume that Zg = 0- Imposing the prior that Ag < 2.2, 
we get the constraint on Ag to be [0.91, 2.20] at 95% 
confidence level. The bounds seem to be sensitive to 
the prior on G. The weak constraint on G is expected 
given the small change in the CMBA power spectra even 
for relatively large variations in G as discussed in Sec- 
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FIG. 6: The marginal distributions of Aq and In Zs obtained 
using the three-year WMAP data set in the instantaneous 
stabilization scenario. Here and thereafter, the maximum of 
the distribution is normalized arbitrarily to 1. 



tion im Now we implement the instantaneous stabiliza- 
tion parametrized by Eq. [51 Since there are already tight 
constraints on the variation of G at redshifts of order 0, 
we impose the prior of In^s > 0. On the other hand, we 
want to constrain the variation of G after CMB decou- 
pling, and thus we impose the prior that Inzg < 6.8. The 
marginal distributions of A§ and In Zg are shown in Fig.[6l 
The 2(7 confidence intervals of Aq and In Zs are [0.95, 1.05] 
and [0, 5.57] respectively. With stabilization, the confi- 
dence interval of Ag shrinks substantially. In Fig. [71 we 
show the contour plot of the joint marginal distribution 
of Aq and In Zg . The triangular shape of the distribution 
is due to the fact that the constraint on Aq is tighter if 
InZs is larger. 

We now turn to the linear stabilization given in Eq. [31 
For smooth variation, we set Zs = 0. The marginal dis- 
tribution of Aq for = is shown in Fig. [51 and the 
resultant 95% confidence interval of Aq is [0.89,1.13]. 

Translating the above constraint on the linear 
parametrization of G to the common form G/G, we get 
G/G = (-9.6 - -h8.1) X 10-12 yr-i. This is comple- 
mentary to the constraints from neutron star mass and 
BBN, which constrain the variation of G in the regimes of 
redshifts ~ 4 and 10^° respectively |8i]. The results are 
summarized in Table [J We see that the CMBA power 
spectra can extend the constraint on the variation of G 
to a large range of redshifts that other experiments and 
observations cannot reach. Improving the low-redshift 
bounds on G helps to tighten the bounds at high red- 



FIG. 7: The contour plot of the joint distribution of Aq and 
In Zs constrained using the three-year WMAP data set in the 
instantaneous stabilization scenario. The inner and outer 
solid lines are the 68% and 95% confidence level contours 
respectively. 
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FIG. 8: The marginal distributions of Aq for the linear sta- 
bilization scenario in Eq. [3]with = and Zs = 0.8 respec- 
tively, constraint using the three-year WMAP data. 



shifts because we can set Zg to a large value. This is 
supported in Fig. [71 for instantaneous stabilization and 
also for linear stabilization in Fig. [51 where we also plot 
the marginal distribution of Aq with Zg = 0.8. The al- 
lowed range of Aq is smaller than that of = 0. 

Although we only consider two types of parametriza- 
tions, they in fact encompass a large class of models in 
which G varies monotonically after photon decoupling. 
If G varies sharply near some redshift z^, this can be 
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TABLE I: The constraints on the variation of G at various 
redshifts. The CMBA constraint fiUs the gap in the "redshift 
ladder" in between neutron star mass and BEN. 





redshift 


G/G yr-i 


Lunar laser ranging [2] 





(1 ±8) X 10"^^ 


Neutron star mass [5] 


0~4 


(-0.6 ±2.0) X 10"^^ 


CMBA (WMAP) 


~ 1000 


(-9.6 ~ 8.1) X 10"^^ 


BBN f6] 


10^° 


(-2.7 ~ 2.1) X 10"" 



approximated by the parametrization Eq. [51 and the 
contour distribution in Fig. [7] can be appUed. On the 
other hand if G varies in a smooth manner, the hnear 
parametrization Eq. [3]can be a good approximation. 

We note that other cosmological parameters that we 
also vary agree well with the standard ACDM cosmolog- 
ical parameters given in [2l|, and so we do not bother 
to write them down. Without stabilization, all cosmo- 
logical parameters are still within 2a from those in the 
WMAP paper. The agreement gets better when we take 
stabilization into account. In these cases, all are within 
Icr. 

From the MCMC runs, we can analyze the degeneracy 
between Aq and other cosmological parameters. We see 
that A§ has some degeneracy with w^, cue dm and rtj. 
The degeneracy with lub and cue dm can be understood 
by the fact that Aq and oji {i = B or CDM) appear in the 
Friedmann equation as the product AgWi. The change in 
relative amplitudes on different scales due to Aq can be 
compensated by changing the relative amplitudes in the 
primordial perturbations characterized by [T3| . As 
mentioned earlier, for instantaneous stabilization, Aq is 
quite strongly degenerate with Inz^. It is well-known 
that the effect of curvature on the spectrum is to shift it 
sideways, somewhat similar to the effect of stabilization 
of G, which is the most important contribution to our 
constraints. Thus we expect that our bounds on Aq will 
be weakened by the inclusion of curvature in the fitting. 

Since the upcoming Planck satellite mission is going 
to probe the temperature power spectrum to as high as 
I ~ 2500 and the E-polarization spectrum to / ~ 1500, 
we expect that there will be tremendous improvement 
in the constraint on Aq. We can forecast the improve- 
ment that Planck will bring us quantitatively using the 
Fisher matrix, which has been widely used to predict 
the expected uncertainties in future experiments (see e.g. 
[isl [22|). Under the assumption of Gaussian perturba- 
tions and Gaussian noise, the Fisher matrix takes the 
form 

where pi is the ith free parameter and Cxi is the Zth 
multipole of the observed spectrum of type X, which 



can be the temperature, temperature-polarization and 
E-polarization spectra. The experimental precision is en- 
coded in the covariant matrix CavixY- We find that with 
the temperature power spectrum alone the current con- 
straint is improved by a factor of 8; when the polarization 
spectrum is included, the bounds will be tightened by a 
factor of 11 relatively to our current bounds. The CMBA 
constraints on G/G will be potentially one of the best 
constraints. Furthermore, it is possible to strengthen the 
constraint by including the matter power spectrum as 
well since enhanced ISW effects are induced by the grav- 
itational stabilization. 



IV. CONCLUSION 

Previously, CMBA was used to constrain the variation 
of G without considering stabilization. Not only are the 
resultant limits weak, but also this assumption does not 
respect many tight local constraints. In this work we con- 
sider two simple and generic parametrizations of G, the 
instantaneous stabilization and linear stabilization. Sta- 
bilization causes appreciable sideway shifts in the CMBA 
power spectra, and hence the sensitivity of CMBA to 
the variation of G is enhanced. We use the three year 
WMAP data to constrain the model parameter (s) and 
other cosmological parameters. For the case of instanta- 
neous stabilization, we simultaneously constrain A§ and 
hvzs to [0.89, 1.13] and [0, 5.57] to 2a intervals. For the 
linear stabilization scenario, Zs is set to and we get the 
95% confidence interval [0.89, 1.13], which is equivalent 
to G/G = (-9.6 - 8.1) X 10~i2 yr-i. Although we only 
concentrate on two types of parametrization, our results 
can be applied to a large class of models in which G varies 
monotonically after CMB decoupling because the varia- 
tion of G in many of these models can be approximated 
by either a step function or a linear function. The con- 
straint derived from CMBA extends the bounds on G up 
to the redshift of about 1000, and so it is complementary 
to other experiments and observations. In particular, it 
fills the "redshift gap" between BBN and the neutron star 
mass constraints. With the forthcoming Planck data, the 
constraints will be improved by a factor of 10 or so, and 
the constraints on G/G from CMBA may be one of the 
best ones. 
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